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Abstract 

We consider the problem of bootstrap percolation on a three dimensional lattice and we 
study its finite size scaling behavior. Bootstrap percolation is an example of Cellular Automata 
defined on the d-dimensional lattice {1, 2, L} d in which each site can be empty or occupied by 
a single particle; in the starting configuration each site is occupied with probability p, occupied 
sites remain occupied for ever, while empty sites are occupied by a particle if at least £ among 
their 2c? nearest neighbor sites are occupied. When d is fixed, the most interesting case is the one 
£ = d: this is a sort of threshold, in the sense that the critical probability p c for the dynamics 
on the infinite lattice Z d switches from zero to one when this limit is crossed. Finite size effects 
in the three-dimensional case are already known in the cases £ < 2: in this paper we discuss 
the case I = 3 and we show that the finite size scaling function for this problem is of the form 
f(L) = const/ In In L. We prove a conjecture proposed by A. CD. van Enter. 
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Cellular Automata are dynamical systems denned on the d-dimensional lattice Z, d in which 
each site of the lattice is occupied by one of finitely many types at each time t. An updating 
rule is defined, which is homogeneous, all the sites follow the same rule, and local, transitions are 



determined by the configuration of types on a finite set of neighboring sites [18, 13 1. 

These models can be thought as interacting particle systems and their connections with statis- 
tical mechanics models have been widely studied in past years (see, for instance, |, 0, |l|, || ) . A 
particular example of cellular automata, known as bootstrap percolation, has been introduced in || 
to model some magnetic systems. More informations on the physical relevance of this model are 
given in 1 10 , |3| . 



In bootstrap percolation only two different types are associated to each site: each site can be 
occupied or not by a particle. In the starting configuration each site is independently occupied with 
probability p, occupied sites remain occupied for ever, while empty sites are occupied by a particle 
if at least £ among their 2d nearest neighbor sites are occupied. The object of primary interest is 
the probability Pf u n(£) that at the end of the dynamics, that is in the infinite time configuration, all 
the sites are occupied. The basic question that has been addressed in physics literature is whether 
by changing the value of the parameter p the system exhibits a sort of phase transition, that is 
whether there exists a critical value p c (£) £ [0, 1] such that if p > p c (£) then p{ u \\{£) = 1, otherwise 
Pfuli(^) < 1. 

Fixed the dimension d, the smaller I the easier empty sites are occupied, hence one expects that 
p c (£) is an increasing function of I. The first rigorous result on this topic is due to van Enter, who 
used the idea of the Straley's argument |1C] to prove that in the case d = 2 and 1 = 2 the critical 



probability is equal to zero. In [[L4|, [T(| Schonmann has proved that p c {£) £ {0, 1}, more precisely 
p c (£) = if I < d, otherwise p c (£) = 1; these results suggest that the most peculiar case is £ = d. 

Before these rigorous results the phase transition scenario in bootstrap percolation models was 
not clear. The technique that had been used to measure the critical probability p c (£) was the finite 



size scaling: a finite volume estimate of the critical probability was found by means of Monte Carlo 
simulations on a finite lattice Al = {1,2, ...L} d , for instance the probability p°jf that one half of 
the samples were completely filled at the end of the dynamics, and the critical value p c (£) was 
extrapolated by means of a suitable scaling function f(L). That is the expression 

p L 5 -Pc(t) L ~°°f(L) (1) 

was supposed to be valid and Monte Carlo data were properly fitted by means of the function f(L) 
(see M and references therein). 

It is rather clear that the estimate of p c strongly depends on the choice of the scaling function 
f(L): the typical choice, when critical effects in second order phase transition are studied, is 
f(L) = const x L~ x l v where v is a suitable exponent. Actually this choice with 1/u = d is correct 
in the case £ = 1, while estimations of p c (£) in the cases 1 = 1 and d > 1 obtained by means of 
Monte Carlo data analyzed through this function f(L) did not fit in the rigorous scenario depicted 
by Schonmann's results [Q]: the problem is that the power law L^ 1 ^ approaches zero too quickly 
and must be replaced by a slower function f(L) = const x (lnL) - ^ -1 ) as suggested by the finite 



volume Aizenman and Lebowitz's results 11]. Indeed the analysis of old and new data performed 
via the correct scaling function yields the correct estimate of the critical probability 0, |^, || . 

In bootstrap percolation on finite lattices is considered in the case £ = 2 and d > 2 and 
it is observed that if p is kept fixed, then in the limit L — > oo the probability to fill {1, L} d 
tends to one whatever the value of p is. But if p — > together with L — > oo, then it is possible 
to find a particular regime in which the probability to fill everything tends to zero. Indeed they 
prove that there exist two constants c + > c_ > such that if p > c^/{\nL) d ~ 1 then pi% — > 1 when 
p — > and L — * oo, while if p < c_/(lnL) d_1 then in the same limit pfJZ — * 0. 

Let us focus on the case d = 3: the choice I = 2 is not the most delicate one, indeed according 
to |l(| even in the case I = 3, that is even in a situation in which it is more difficult to fill empty 
sites, the critical probability is still zero. Hence one can guess that in the three-dimensional case 



if I = 3 then the finite scaling function is no more the Aizenman-Lebowitz one, but a function 
approaching zero more slowly. Our aim is to study this case and to show that results similar to 
those in H, and conjectured by A. CD. van Enter, can be proved if the scaling function is replaced 



by /(-£/) = const/ lnlnL. This problem has been proposed in 16] as Problem 3.1; we notice also 
that related problems have been discussed in [^, [l5|] . 

An interesting follow up would be the generalization of our results to the d-dimensional case by 
considering arbitrary dimension d and £ = d. In this case one expects f(L) = const/ (In In ... InL) 
with the logarithm applied d — 1 times. 

We next define the particular model of bootstrap percolation that we are going to study and we 
introduce some notations. Let us consider the lattice Z 3 and the discrete time variable t = 0,1, 2, .... 
To each site x G 71? we associate at each instant of time t a random variable Xt(x) which takes 
values in {0, 1}; depending on X t (x) equals or 1 we say that the site x is empty or occupied. 
We denote by O := {0, l} z3 the space of configurations and by Xt G £1 the configuration of the 
system at time t. The initial configuration Xq is chosen by occupying independently each site of 
the lattice with probability p (initial density). Then the system evolves according to the following 
deterministic rules: 

• if X t (x) = 1, then X t +i(x) = 1 (l's are stable); 

• if X t (x) = and x has at least three occupied sites among its six nearest neighbors, then 
X t+ i(x) = 1; 

• Xt + i(x) = otherwise. 

We omit in our notations the dependence of the process Xt on the initial density p; when the initial 
density will be different from p it will be clearly stated. 

The primary object of interest for this problem is the final configuration 



X := lim X t , (2) 

t^oo 



that is the configuration attained by the system when the dynamics stops. 

Let us consider a subset A C Z 3 , we denote by X\ it the process restricted to A with free 
boundary conditions, that is without taking into account sites outside A. In this case Xa,*, X\ ; t(x) 
and X\ will respectively denote the configuration at time t, the value of the random variable at 
time t and site x, and the final configuration. 

Given two arbitrary sets Ai,A2 C Z 3 we will consider the process X^ t restricted to Ai and 
with occupied sites in A2 (the sites in Z 3 \ (Ai U A2) are not taken into account). We will obviously 
omit Ai in the notation if Ai = Z 3 , while we will omit A2 if A2 = 0. 

Definition 1 Following ffij we say that a set A C Z 3 is internally spanned if it is entirely covered 
in the final configuration of the dynamics restricted to A, that is if 

G A X A (x) = 1 . (3) 

We have now introduced the bootstrap percolation model in dimension d = 3 and with pa- 
rameter I = 3; as we have already stated, this is the most delicate three-dimensional bootstrap 
percolation model. Indeed £ = 3 is the highest value of the parameter i such that the critical prob- 
ability in infinite volume is equal to zero, that is for each positive initial density p the probability 
that the whole lattice will be completely occupied by particles at the end of the dynamics is equal 
to one. 

In this particular case we examine the question of finite size scaling, that is, following Q, we 
consider the process X\ L j on the finite cube Al = {1, ...,L} 3 of size L and we perform the limit 
L — > 00 and p — ► 0. We prove that there exists a particular regime in which the probability that 
the cube Al is internally spanned is zero in the limit L — > 00 and p — > 0. 

Theorem 2 Let us consider the cube Al of side length L and the process X% with initial density 
p; let us denote by R(L,p) the probability that the cube Al is internally spanned. There exist two 
constants c+ > c_ > such that 



• R(L,p) — ► 1 if (L,p) — > (oo,0) in the regime p > c + /lnlnL 

• R(L,p) — ► if (-L, p) — ► (00, 0) in the regime p < c_ / In In L 

We start to prove the first part of Theorem [2f this is the easy part of the theorem and its 
proof has already been sketched in jjj]. The idea of the proof relies on the notion of critical length: 
suppose p is small, if an occupied cube has size large enough, then the probability to find on its 
faces two-dimensional occupied square droplets so large to cover the whole faces of the cube (two- 
dimensional critical droplets) is close to one. Obviously if p goes to 0, then the size of the cube 
must diverge and this must happen fastly enough; in this way one will find that the critical length 
is of order exp (const /p). Then, one estimates that in order to have a high probability of finding a 
critical droplet, the size L of the cube must be of order exp exp (const /p). 

Now, we come to the proof: from results in [16] one easily obtains that there exists a constant 



c\ > such that, given a cube A;, if p > c\/\nl then there exists a constant a\ > such that 

00 

P(A ; covers 1? \ A t occupied at t = 0) > JJ(1 - e" aifc ) . (4) 

k=l 

We consider a large cube Al and we estimate R(L,p): 

R(L,p) > P(m Al 3Ai occupied at t = 0, A/ covers Z 3 ) = 

P( Ai covers 1? \ in Al 3A^ occupied at t = 0) P(in Al 3A^ occupied at t = 0) (5) 

Now, from equation (^j it follows that if I > exp(ci/p) then the first factor in (||) tends to 1 when 
p goes to 0. Then, we just have to prove that with such an I it is possible to choose L such that 
the second factor of @ tends to 1, as well. Indeed by partitioning the cube Al in disjoint cubes of 
size I one has 

P(in Al 3Ai occupied at t = 0) > 1 - (1 - j/) L/z ; (6) 

by choosing c + > 3ci > and L > expexp(c + /p), the right hand term in the equation above tends 
to 1 when p goes to 0. This completes the proof of the first part of Theorem 0. 



Next we prove the second part of Theorem ^ by finding a suitable upper bound to the probability 
that a cube Ax, is internally spanned. First of all we give some more definitions: given a site i£Z 3 
we denote by (x%, x%, X3) its three coordinates and given a set A C Z 3 we define its diameter 

d(A) := sup{\ Xi -yi\ : x G A, y £ A, i £ {1,2,3}} , (7) 

that is d(A) is the side length of the minimal cube surrounding the set A. We say that A C Z 3 is 
a region of Z 3 if and only if it is nearest neighbors connected. We note that if Ai and A2 are two 
regions of Z 3 and Ai U A2 is a region as well, then d{A\ U A2) < e?(Ai) + &{k-i). 

Now we adapt to our situation a key Lemma of || that describes what happens at a smaller 
scale inside an internally spanned region. 

Lemma 3 Let Ai be a region of 71? . If Ai is internally spanned, then for all k such that 1 < k 
and 2k + 1 < d(A\) there exists at least a region A2 included in A\ which is internally spanned and 
such that k < d(A2) < 2k + 1. 

Proof. We build the configuration X\ t by means of the following algorithmic procedure. Let Co 
be the collection of the sites occupied at time zero. Suppose we have built a collection of regions 
internally spanned C n , we define a rule to build the collection C n+ \: 

• If there exist two regions A, B of C n such that A U B is still a region, then we set 

C n+1 :=C n U{AiJB}\{A,B} , (8) 
that is is obtained by replacing in C n the two elements A and B by A U B. 

• If no such regions exist then we choose a site x not belonging to any set in C n and having three 

neighbors in the set (J A. We denote by Ai with 1 < i < r the r regions of C n containing 

Aec n 

a neighbor of x, and we set 

C n+ i :=C„U j|J AiU{x}| \{A u ...,A r } . (9) 

7 



• If no such site x exists, the algorithm stops. 

Notice that for each n, the regions of C n are internally spanned. Since Ai is internally spanned the 
procedure ends for some m such that C m = {Ai}. Moreover we have that max{d(A) : A G Co} = 1, 
max{d(A) : A G C m } = d(Ai) and for any n < m — 1 

max{d(A) : A G C n+1 } < 2max{d(^) : A G C n } + 1 . (10) 

Hence there exists n such that: 

k < max{d(A) : A G C n } < 2k + 1 , (11) 

which means that in C n there is an internally spanned region A of diameter d(A) such that k < 
d{A) < 2k + 1. ■ 

Definition 4 Lei us consider a cube A in Z 3 . M^e say t/iat A is crossed, or that there is a crossing 
in A, if and only if in the final configuration X\ of the dynamics restricted to A there is an occupied 
region joining two opposite faces of the cube A. 

We note that for any region A the following inclusion stands: 

{^4 is internally spanned} C {the smallest cube surrounding A is crossed} . (12) 

Hence, for any L and any k such that 2k + 1 < L, we have: 

R(L,p) < P(3l, k < I < 2k + 1, 3A Z C A L , A/ is crossed) < 



(k + 1) L 3 max P(A ; is crossed) (13) 

k<1<2k+1 



Thus, 



R(L,p) < L min (k + 1) max P(A t is crossed) . (14) 

1<k<(L-1)/2 k<K2k+1 



We have reduced the estimate of R(L,p) to the estimate of the probability that a cube A; is 
crossed and by symmetry we can consider the case of a crossing along the first reticular direction 
(denoted by e\ in the sequel): 

P(Ai is crossed) < 3 P(A; is crossed along ei) . (15) 

In order to estimate the second hand term of equation (|l^) we reconduct to a two-dimensional 
situation by properly cutting the cube A; in slices of thickness two and perpendicular to the first 
reticular direction. For the sake of definiteness we suppose A^ := {1, 2, I} 3 , I an even number 
and we define the slices 

Tk := {i £ A| : x\ = 2k — \ or x\ = 2k} , 1 < k < l - . (16) 

We define a map s associating to each site in a slice the only nearest neighbor along the first 
reticular direction belonging to the same slice: 

VxGA, ,(*):=( & + ^ *i ^ odd 

[ \x\ — 1, xii X3) if x\ is even 

The process X T t k , restricted to the slice Tjt and with all the sites in A; \ occupied, dominates 
the original process in the same slice: 

Vk G {1, h Vx G T k Vt > X t {x) < x£ A t Tfe (x) . (18) 
2 ^ ' 

In each slice for k in {1, 1/2} we define a new process Y t k . 

Definition 5 We consider all the sites in Ai \ T k occupied and define the process Y t k on {0, l} Tfe 
as follows: for any x in 

. Y k (x) = max(X^ J k (x) , X% J k (s(x)))= max(X (x) , X (s (x) ) ) 

• ifY t k {x) = 1 then Y k +l {x) = 1 

• if'Y k (x) = and x has at least 3 occupied sites among its 6 nearest neighbors, then Y k +1 (x) = 1 
and Y t k +1 (s(x)) = 1 

n 



• Y k +l (x) = otherwise 

The mechanism to build Y k +1 (x) is the one used for followed by an additional step increasing 

the configuration. This mechanism ensures that for all t and any x in one has Y k (x) = Y k (s(x)). 
Finally we introduce a two-dimensional process that will be used to estimate fll5|). 

Definition 6 Let us associate to each slice for k in {1, ...,1/2} a two-dimensional I x I square 
Q k := {1, 2, I} 2 ; then on each square Q k we define a process Zqk t by 

V(x 2 ,x 3 )eQf Z Qf , t (x 2 ,x 3 ) :=Y k (2k-l,x 2 ,x 3 ) = Y k (2k,x 2 ,x 3 ) . (19) 

The processes Zqk t , 1 < k < 1/2, are independent and they are two-dimensional bootstrap perco- 
lation processes with parameter 1 = 2 and initial density q = 2p — p 2 . Furthermore these processes 
dominate the original process on the slices in the sense: 

Vk £ {1, \-} Vx = (xi,X2,xs) e T k Vi > X Tk ,t{x) < Z Qk t (x 2 ,x 3 ) . (20) 

2 ' ' 

The two-dimensional processes can be used to estimate the probability that a cubic region is 
crossed: we consider the (1/2) x I x I parallelepiped Vi obtained by collecting the 1/2 squares Qf 
and we denote by Z-p t the configuration on V[ defined as follows: 

Vxi G {1, l -} Vx2,x 3 G {1, ...,/} Z Vl (xi,x 2 ,x 3 ) := Zq*i(x2,x 3 ) , (21) 

where Zqk with k in {1, ...,1/2} is the final configuration of the process Zqk t . Finally we have the 
way to estimate ([l5[): 

P(Ai is crossed along e±) < P(in Zp l there is a crossing along e{) . (22) 

Now we consider the two-dimensional bootstrap percolation with parameter 1 = 2 and initial 
density q. We define as well the concept of "being internally spanned" and we denote by S(l) a 
square of side length I. We recall that the final configuration, for such a process, is a union of 
separated rectangular regions and we state a few results: 

1 n 



Lemma 7 (Aizenman - Lebowitz []3|]) For all k > 1, a necessary condition for S(l) to be in- 
ternally spanned, where k < I, is that it contains at least one rectangular region whose maximal 
side length is in the interval [k, 2k + 1] which is also internally spanned. 

Proof of Lemma [7] is given in Q . 

Lemma 8 Let A be a rectangular region of side lengths l\ and I2; suppose l\ < l 2 - For q small 
enough one has 

P(A is internally spanned) < (4/ 2 <?)' 2/2 . (23) 

Proof. If A is partitioned in Z2/2 disjoint slabs of width 2, a necessary condition for A to be 
internally spanned is that each slab contains initially an occupied site. Hence: 

P(A is internally spanned) < ^1 — (1 — q) 2ll ^j ^ < 

exp (| In (1 - exp(2/ 2 ln(l - q)))j < 

exp(! In (-2/ 2 Hl-q))\ (24) 
For q small enough, ln(l — q) > —2q, whence 

P(A is internally spanned) < exp f -| \n{<U 2 q) J = {Al 2 q) h/2 . ■ (25) 



Lemma 9 For any I and any k < (I — l)/2 let £ be the event: S(l) contains a rectangular region 
internally spanned whose maximal side length belongs to the interval [n, 2k + 1]. For q small enough 
one has 

P{£) < I 2 (2k + l) 2 exp exp(-4(2K + l)qfj . (26) 



Proof. We suppose q small enough to have ln(l — q) > —2q and we bound the probability of 
the event £ as follows. 



< / 2 (2k + l) 2 max exp ^ ln(l - exp(2/ 2 ln(l - q))) < 



V (2k + l) 2 exp ^- ln(l - exp(2(2K + 1) ln(l - g)))J < 
/ 2 (2k + l) 2 exp Q ln(l - exp(-4(2K + < 

Z 2 (2k + 1) 2 exp(-^ exp(-4(2« + . ■ (27) 

Now we come back to the proof of the upper bound. Let us consider a > 0, we notice that 
P < q < 2p. We denote by 1 a configuration in a square Qf with all the sites occupied. We still 
increase the configuration Zgk by setting Zgk = 1 in case Zqk contains at least one rectangular 
region of maximal side length larger than a/q. Supposing I > a/q and q small enough so that 
a/q > 3, by applying Lemma ^ with k = a/3q one has 

^-D^ (| + ^ e*p (-£ exp (-4 (| + 
l 2 a 2 ( a . M \ 

~T exp v% exp( ~ 4a) J (28) 

We suppose that a is small enough to have exp(— 4a) > 1/2; thus, 

^ Qf =I)<^exp(--) . (29) 

Let M be the number of indices k such that one has Zgk. = 1 and let k(l), ...,k(M) be these 
indices arranged in increasing order. Let £\ be the event {there is a crossing along e\ in Z-p^. We 
decompose this event as follows: 

1/2 

P(£ 1 ) = P(S l ,M = 0) + J2 P(£u M = m,k(l) = h,...., k(m) = i m ) . (30) 

m=l ii<-"<i TO 

Let i < j be two indices in {1, ...,l/2}; by £(i,j) we denote the following event: there exists a 
sequence of disjoint rectangular regions (Rh, 1 < h < H) in Z-p l such that 

i o 



• Ri is included in Q\, Rh is included in Qj; 

• the regions Rh with 2 < h < H — 1 are included in |J Q^; 

i<h<j 

• the maximal side length of all these regions is strictly less than a/q: we denote by r/j the 
maximal side length of R^; 

• for each h, 1 < h < H — 1, a site of Rh is the neighbor of a site of Rh+i, 

• all the sites of these regions are occupied in Z-p t . 

We remark that H is free, however it has to be larger than j — i + 1. Moreover the sequence of 
rectangles (Rh, 2 < h < H — 1) can go back and forth between the squares Q\ +l and Qj" 1 - We 
make the convention that for any i the events £(i,i) and £(i,i — 1) are the full events. We have 
the following estimate: 

P(£i, M = m, k(l) = i±, k(m) = i m ) < P(M = m, k(l) = i±, k(m) = i m , 

6(1, ii - 1), 6(h + l,i 2 - I),--, 6(i m -i + l,i m - 1), 6 (i m + 1, l -)) = 

m+l 

n p(6(i h _ 1 +i,i h - 1)) p(z Qh =!)••• p(z Q]m = i) < 
h=i 

exp v^q)) n + - x )) ( 3i ) 

where we have set i$ = and ? m +i = | + 1. 

In order to estimate P(6(i,j)) with i < j, we consider a fixed sequence n, r# in {1, a/g — 1}; 
the number of sequences R±, ...,Rh of rectangles with maximal sides r\, ■■■,ru and satisfying the 
above requirements is smaller than 

l 2 r x x r\r 2 x r|r 3 x ... x r\_ x r H < I 2 (r x r 2 ■ ■ ■ r H f . (32) 

Notice that several rectangles of the sequence R\, Rh can belong to the same slice. However the 
rectangles are disjoint and the events that they are internally spanned depend only on the dynamics 

1 Q 



restricted to the rectangles, hence these events are independent, so that 

P(R\, ■■■,Rh are occupied in Zp t ) < P(Ri, ■■■jRh internally spanned) 

P(R\ internally spanned) • • • P(Rh internally spanned) < (4rig) ri//2 • • • (4r# 

where in the last inequality we have used the Lemma ||. Thus 

P(£(i,j))= Yl P£Ri,-,Rh realizing £(i, j)) < 
H>j-i+l 

E l 2 (ri'- r H f (4r iq y^ ■ ■ ■ (4r H q) r »/ 2 = 
H>j—i+l n,...,rH<a/q 

H 

£ I 2 ( £ r^Arq)^ 

H>j-i+l \l<r<a/q 

We estimate the sum ^ r 3 (4rq) r / 2 as follows: 

l<r<a/q 

r 3 {Arqfl 2 = ^ r 3 {Arq) r ' 2 + ^ r 3 (4rg) r / 2 < 

l<r<a/ij l<r<8 9<r<a/g 

\ 4 

a 



8 3 (32g) 1 / 2 + - max {Arq) r ' 2 . 

\qj 9<r<a/q y 

Let f(r) = (4rq) r ^ 2 ; for a small f(r) is decreasing on [9,a/q[, whence 



£ r 3 (4rg) r / 2 < 8 3 (32g) 1 /2 + a 4 36 9 /y/ 2 < b Q q^ 2 
l<r<a/q 

where bo is a constant not depending on a. Thus, 



P(£(i,j))< £ z 2 (W^) H 

Finally, for q small enough, so that bo^/q < 1/2, one has 



P(S(i,j)) < 2l 2 (boVqy- i+1 = 2/ 2 exp((j -i + 1)(Ih6q + ^lng)) 



Coming back to inequality (31) 



P(£i, M = m, k(l) = h, fc(m) = i m ) < 



Pa 2 

2~ eX P 



a 



(2l 2 ) m+1 exp 



and 



12qJ 



1/2 



I \ ( 1 

m] ( In bo H — In < 



(39) 



P(£ 1 ) = P(£ 1 ,M = 0) + Y J Yl P(Sx, M = m,k(l) = h,...., k(m)=i m )< 

1/2 



m=l ii<— <i m 



(2/ 



2/ 2 exp 



A(ln6 + ilng 



' /2 W 



£ 

m=0 



exp 



i-m) ( ln fro + ^ In < 



< 



-^expi-— ,ex P 



ln bo H — ln 



(40) 



This is our estimate of the probability of a crossing in the case I > a/g. On the other hand, in the 
case I < a/q we estimate the probability of a crossing by considering directly the event £(1, 1/2): 



P(Ai is crossed along e\) < P[£(l, -)) < 2l 2 exp 



I ( 1 
- In6 + -ln, 



(41) 



We note that it would have not been possible to use the same strategy in the case I > a/q because 
the probability of having a very large internally spanned rectangle does not vanish. 
Supposing that I < exp(a/120g) one has 



l 5 a 2 

2~ eX P 



a \ 



ln b + - In q 



-^ eXP {-2Tq- lnb °- 1 2 lnq 



(42) 



Now, if q is sufficiently small so that the right hand term is smaller than 1 and ln bo < —(1/4) ln q, 
under the hypothesis 2 < I < exp(a/120g), from @, @ and © one has 



P(Ai is crossed) < 6l 2 exp 



l + a < 6 j3 exp n lnq 



(43) 



Hence, there exists a > such that for p sufficiently small and 2 < I < exp(a/240p) 



P(Ai is crossed) < 6/ a exp ( - ln(2p) 



(44) 



Finally, we can use equations (14) and ( |44| ) to estimate the probability R(L,p) and to complete 
the proof of Theorem [2|. First we consider the case L < exp(ao/p) with ao = a/240 and we write 



R(L,p)<L 3 min (k + 1) x 6(2k + l) 3 exp ( J ln(2p)) < 6L 7 exp ( ^- ln(2p) ) 
i<n<(L—i)/2 \8 / y24 J 



(45) 



and we remark that the right hand term goes to zero in the limit p — > and L — > oo. On the other 
hand, in the case L > exp(ao/p) one can restrict the minimum to 2k + 1 < exp(«o/p) and write 



R(L,p) < 6L 3 exp 



( 4ap 



cxp 



1 



2 1 ' \-) H2P) 



(46) 



From the estimate above it is clear that there exists a positive constant c_ such that if L is less 
than expexp(c_/p) then R(L,p) goes to in the limit L — > oo and p — > 0. This completes the 
proof of Theorem S. I 



Acknowledgements 

The authors express their thanks to A. CD. van Enter who suggested the topic of this paper and 
to R.H. Schonmann for useful suggestions and comments. 



1 « 



References 

[1] J. Adler, A. Aharony (1988). Diffusion percolation: I. Infinite time limit and bootstrap perco- 
lation. J. Phys. A: Math. Gen. 21, 1387. 

[2] J. Adler, D. Stauffer, A. Aharony (1989). Comparison of bootstrap percolation models. J. 
Phys. A: Math. Gen. 22, L297. 

[3] M. Aizenman, J.L. Lebowitz (1998). Metastability effects in bootstrap percolation. J. Phys. 
A: Math. Gen. 21, 3801. 

[4] N.S. Branco, R.R. Dos Santos, S.L.A. de Queiroz, J. Phys. C 17, L373 (1984); M.A. Khan, 
H. Gould, J. Chalupa, J. Phys. C 18, L233 (1985); N.S. Branco, S.L.A. de Queiroz, R.R. Dos 
Santos, J. Phys. C 19, 1909 (1986). 

[5] J. Chalupa, P.L. Leath, G.R. Reich (1979). Bootstrap percolation on a Bethe lattice. J. Phys. 
C: Solid State Physics 12, L31. 

[6] A. CD. van Enter (1987). Proof of Straley's argument for Bootstrap Percolation. J. Statist. 
Phys. 48, 943. 

[7] A. CD. van Enter, J. Adler, J. A. M.S. Duarte (1990). Finite-size effects for some bootstrap 
percolation models. J. Statist. Phys. 60, 323. 

[8] A. CD. van Enter, J. Adler, J. A. M.S. Duarte (1991). Addendum: Finite-size effects for some 
bootstrap percolation models. J. Statist. Phys. 62, 505. 

[9] D. Griffeath (1998). Cyclic random competition: a case history in experimental mathematics. 
Notices of the American Mathematical Society 35, 1472. 

[10] P.M. Kogut, P.L. Leath (1981). J. Phys. C: Solid State Physics 14, 3187. 



[11] R. Le Normand, C. Zarcone (1984). "Kinetics of Aggregation and Gelation", eds. F. Family 
and D.P. Landau (Elsevier, Amsterdam). 

[12] T.S. Mountford (1995). Critical length for semi-oriented bootstrap percolation. Stochastic 
Process. Appl. 56, 185-205. 

[13] J. von Neumann (1966). Theory of self-reproducing automata. Univ. of Illinois Press, Urbana. 

[14] R.H. Schonmann (1990). Critical Points of two-dimensional bootstrap percolation-like cellular 
automata. J. Statist. Phys. 58, 1239. 

[15] R.H. Schonmann (1990). Finite size scaling behavior of a biased majority rule cellular automa- 
ton. Phys. A 167, 619-627. 

[16] R.H. Schonmann (1992). On the behavior of some cellular automata related to bootstrap 
percolation. Ann. Probab. 20, 174. 

[17] T. Toffoli, N. Margolus (1987). Cellular automata machines. A new environment for modeling. 
MIT Press. 

[18] S. Ulam (1950). Random processes and transformations. Proc. Internat. Congr. Math., 264- 
275. 

[19] G.Y. Vichniac (1984). Simulating physics with cellular automata. Phys. D 10, 96. 

[20] S. Wolfram (1983). Statistical mechanics of cellular automata. Rev. Modern Phys. 55, 601 
(1983). S. Wolfram (1986). Theory and applications of cellular automata. World Scientific, 
Singapore. 



1 8 



